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The discovery of leptonic CP violation is one of the primary goals of next-generation neutrino 
oscillation experiments, which is feasible due to the recent measurement of a relatively large leptonic 
mixing angle ^13. We suggest two new working observables AA™^ = max[A^^(^)] — min[A^^(5)] 
and /\A^{5) = A^{5) — A^J(O) to describe the CP- violating effects in long-baseline and atmo- 
spheric neutrino oscillation experiments. The former signifies the experimental sensitivity to the 
leptonic Dirac CP-violating phase 5 and can be used to optimize the experimental setup, while the 
latter measures the intrinsic leptonic CP violation and can be used to extract 5 directly from the 
experimental observations. Both analytical and numerical analyses are carried out to illustrate their 
main features. It turns out that an intense neutrino beam with sub-GeV energies and a baseline of 
a few 100 km may serve as an optimal experimental setup for probing leptonic CP violation. 
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I. INTRODUCTION 

The phenomenon of neutrino oscillations in vacuum 
and matter can be described by six fundamental parame- 
ters: three lepton flavor mixing angles {^127 ^137 ^23} 5 ^^o 
neutrino mass-squared differences {Am2i, Amg^Ll? 
one Dirac-type CP-violating phase 5. Due to a number 
of elegant neutrino oscillation experiments in the past 
decades, both {^127^^21} {^235 l^^iill have been 
measured with reasonably good accuracy [1 . Until re- 
cently, the smallest leptonic mixing angle ^^^3 has been 
found to be relatively large in the Daya Bay [2 and 
RENO [3 reactor neutrino experiments. This great dis- 
covery enhances the capability of the next-generation ex- 
periments to pin down the neutrino mass hierarchy (i.e., 
the sign of /S.m\i) and eventually to determine the lep- 
tonic Dirac CP-violating phase. 

An important question is how to characterize the lep- 
tonic CP-violating effects in neutrino oscillation exper- 
iments. For neutrino oscillations in vacuum, it is evi- 
dent that the CP asymmetry, usually defined as A^ = 
■^cx(3 ~ ^a(3 ^ ^^^^ for a 7^ /3, is well determined by (5, 
where P^^ = P{v^ v^) and P^^ = P(v^ v^) stand 
for the neutrino and antineutrino transition probabili- 
ties, respectively. For long-baseline neutrino oscillation 
experiments, however, matter effects [H |5] can be sig- 
nificant and induce fake CP violation, since the Earth 
matter itself is CP asymmetric. In this case, the in- 
trinsic CP violation due to 8 in A^ is obscured by ex- 
trinsic CP violation caused by matter effects. On the 
other hand, one can determine 5 by just measuring the 
probabilities P^o as precisely as possible. In this case. 
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if one defines AP^J((5) = P^^{8) - P^^(O), the fake 
CP violation can be removed, since AP^J(^) vanishes 
for 5 = 0. Furthermore, it has been suggested [6 that 
AP^ = max[P^^((5)] — min[P^^((5)] can be utilized to 
quantify the experimental sensitivity to (5, where the 
maximum and minimum are obtained by freely varying S 
in [0, 27r). Both AP^J((5) and AP^^ have been studied in 
detail by using neutrino oscillograms of the Earth [71 18]. 

Since the description of leptonic CP violation should 
reflect the difference between neutrinos and antineu- 
trinos, we suggest AA^P(^) = A^^{6) - ^^^(0) and 
AA^^ = mdix[A^^ (S)] — mm[A^^ (S)] as working observ- 
ables to signify the intrinsic CP violation and the experi- 
mental sensitivity to S. First, we make a general compar- 
ison among; all five quantities (i.e., A^^, AP^J, AP^^, 
AA^^, and AA^^) and point out their advantages and 
disadvantages in describing leptonic CP violation. Then, 
we perform a detailed numerical study of them, and de- 
scribe their main features by using approximate analyt- 
ical formulas. Finally, we investigate the experimental 
setup that is optimal for the determination of 5. 



II. MEASURES OF LEPTONIC CP VIOLATION 

First of all, we briefly review the formulation of three- 
flavor neutrino oscillations in matter, which is rele- 
vant for long-baseline experiments. The flavor tran- 
sition of neutrinos propagating in matter is governed 
by the effective Hamiltonian H^^ = ^ V, where 
= U ■ diag(0, A21, A31) • [/t with A21 = AmlJ2E 
and Ag]^ = /S.m\i/2E being the low and high oscillation 
frequencies and E is the neutrino energy. The matter 
potential is F = >/2Gpngdiag(l, 0, 0), where is the 
electron number density. In the standard parametriza- 
tion, the leptonic mixing matrix is = 50i2X\Oi2i 
where O^^ denotes the rotation in the i-j plane with an 
angle 0^^ (for ij = 12, 13, 23) and F^ = diag(l, 1, e^^) with 
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5 being the leptonic Dirac CP-violating phase. Since V 
is invariant under any rotations in the 2-3 plane, one can 
transform to flavor basis [v^^v^^i)^)^ = ^23(^65 ^/^^ ^r)"^ 

with U22, = 02^^^ such that H'^^ = ^23^eff^23 is i^" 
dependent of ^23 S in this basis. The amplitude of 
neutrino flavor transition is A{i'^ v^) = S^^^ where 
the evolution matrix Six) satisfies the Schrodinger-like 
equation i[d6'(x)/dx] = H^^{x)S{x) with the initial con- 
dition 6'(0) = 1. For constant matter density, we have 
S{x) = exp(— ii^g^x) with x being the distance that neu- 
trinos propagate. If the evolution matrix correspond- 
ing to H^^ is denoted as S\ we have S = ^23 ^'^23- 
The oscillation probabilities of neutrinos are given by 
P^l^ = \A{h'^ ^/3)P = I^/3q;P' while those of antineu- 
trinos P^^ can be derived from the same effective Hamil- 
tonian but with the replacements S ^ —S and V ^ —V. 

As shown for example in Ref . [9 , two out of nine prob- 
abilities P^^ are independent. Now, we choose them as 
P^g and P^^. Our choice of P^g and P^^ has two advan- 
tages: (1) Most of the neutrino detectors are designed 
for detection of electrons and muons and their antipar- 
ticles; (2) Both appearance and disappearance channels 
are included, which have very different sensitivities to the 
neutrino mass hierarchy and 5. 

We further show the dependence of P^g and P^^ on S. 

Since S = f^23'S''^23 5 where S' is independent of ^23 ^ind 
^, it is straightforward to prove that P^g and P^^ can be 
written as [10^ 

P^g = a cos ^ + 6 sin ^ + c , 

Pfifi = /cos(5 + ^cos2(5 + /i , (1) 

where the relevant coefficients {a, 6, c} and {f^g^h} are 
independent of 6 and their exact expressions can be found 
in Ref. [11]. Using the two independent probabilities in 
Eq. ([1]), 

one can readily find the exact expressions for all 
the other probabilities. 

For neutrino oscillations in vacuum, the measure of 
leptonic CP violation can be taken to be the Jarlskog in- 
variant J = 5;l2<^12'^23^23'^13<^13 ^5 whcrC S^j = SluO^j 

and a-. = cos6>,.. In fact, the difference between neu- 
trino and antineutrino probabilities is P^^ — P^,^ ex J' for 
a ^ p. For neutrino oscillations in matter, however, the 
situation is complicated by the CP-asymmetric medium, 
since only particles rather than antiparticles are present 
in Earth matter. In this case, the discrepancy between 
Pq,^ and Pq,^ receives contributions both from S and mat- 
ter effects. Therefore, it is natural to ask which measure 
is the best to extract information on 6 from observations, 
and to find the optimal experimental setup to measure S 
in future neutrino experiments. In the literature, there 
already exist three distinct measures: 

(i) ^S^(^) = ^afdi^) ~ ^afsi^) denotes the differences 
between neutrino probabilities Pq,^(^) and antineutrino 
probabilities P^^{d)^ where the dependence on 6 is explic- 
itly displayed. In long-baseline experiments, where mat- 
ter effects play an important role, A^^{5) is no longer the 



best measure of intrinsic CP violation, since ^ 

even for (5 = due to matter effects. With help of the ex- 
act formula for P^g in Eq. ([T]) as well as the counterpart 
P^g with coefficients {a, 6, c} in the antineutrino channel, 
we obtain 

^M^(^) = Aacos(5 + A6sin(5 + Ac , (2) 

where Aa = a — a; likewise for Ab and Ac. Obviously, 
A^f{S) follows the same dependence on S as P^g((^) and 
P^^{5). Note that we focus on the appearance chan- 
nel ^ei^e)^ disappearance channel 
^fii^ij.) ~^ ^fii^fi) can be discussed in a similar way. 

(ii) AP^^{5) = P^^{5) - P^^(O) denotes the differ- 
ences between the probabilities P^^ for an arbitrary 5 
and those for (5 = 0. This measure is intended to re- 
move the fake CP violation induced by matter effects, 
which has the advantage that only the neutrino channel 
is involved. Using Eq. (IT]), we find 



AP^^{5) = 2^a^^h^ 




with tancj = h/a. It is now evident that AP^^ (S) is 
proportional to sin((5/2), and vanishes for 5 = 0. 

(iii) AP^^ = max[P^^((5)] - mm[P^^{S)] denotes the 
variation of the probabilities P(^p{S) for 6 varying in 
[0, 27r). Such a measure is actually independent from the 
true value of 5, which is yet unknown, so it should be 
useful to find an optimal experimental setup that is most 
sensitive to ^. In fact, we have 

AP;^g = 2Va^^b^ . (4) 

Hence, one can observe that AP^ essentially determines 
the amplitude of AP^^{S). The basic strategy to probe 
S may first be to optimize the experimental setup with 
help of AP^ , and then to extract 6 from the observation 
of AP^gP((5). In this sense, AP^g and AP^^{S) can be 
grouped as a pair of working observables to probe S. Both 
of them are based on P^g. Similarly, one can consider 
AP^ and AP^^ (S) in the antineutrino channel. 

We argue that the proper measures of leptonic CP vi- 
olation should manifest the difference between neutrinos 
and antineutrinos. In principle, most of the proposed 
long-baseline experiments are equally operative in the 
neutrino and antineutrino channels. Therefore, based on 
Eq. we suggest a new pair of working observables: 

(iv) AA^J(J) = A^^ (S) - A^^ (0) signifies the intrinsic 
CP violation, compared to A similar quantity, 
but with a different normalization, was previously con- 
sidered Jl2, 13 . Using the exact expression for A^f (S) 
in Eq. p|, we arrive at 

A^^P(^) = 2V(Aa)2 + (A6)2sin ^ sin (^oj' - (5) 

with tancj' = Ah/ Aa. Note that this result takes the 
same form as that of AP^f {5) in Eq. ([3|, except for the 
relevant coefficients. 
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(v) AA^^ = max[A^P((5)] - mm[A'^^{6)] denotes the 
variation of A^^{6) for 6 varying in [0, 27r). The extrinsic 
CP-violating effects cancel in AA^^. It is straightforward 
to show that 



AA- = 2v/(Aa)2 + (A5)2 . 



(6) 



The size of AA^^ determines the magnitude of AA^f{S) 
through the ^-independent coefficient. 

We expect that AA^^ and AA^^ can be implemented 
to extract and to optimize the experimental setup, sim- 
ilar to AP^J and AP^. Nevertheless, the former ones 
contain the difference between neutrino and antineutrino 
probabilities, so these two sets of measures are not equiv- 
alent. In the following, we will present approximate and 
analytical results for all the above measures of leptonic 
CP violation, and the numerical results as well. Further- 
more, the optimal experimental setup for probing S is 
considered and compared with the ongoing and upcom- 
ing neutrino oscillation experiments. 



III. ANALYTICAL & NUMERICAL RESULTS 

We define a A^i/A^-^, A A^-^L/2 with L be- 
ing the distance between the source and detector, and 
A = V/A^i^ where a denotes the ratio of the low and 
high oscillation frequencies, A is the phase correspond- 
ing to the high oscillation frequency, and A stands for the 
strength of matter effects. The probabilities in matter of 
constant density have been calculated in Ref. [9 to the 
second order in both a and Si^. According to the cur- 
rent neutrino oscillation data, we have a ^ V^sf^ ^0.03 
[MHTB]. Therefore, one can expand the probabilities in 
terms of a and 5^3, and neglect all higher-order terms of 
O(a^). Using the approximate formula for P^^, we can 
identify the corresponding coefficients 

^ sinAA sin(A - 1)A 
a ^ +8aJr : r^cosA 



A A-1 
sin AA sin(A - 1)A 



A A-1 

2/ 



sin A 



4,2 2 sin'{A-l)A 



(7) 



^3 



S12C12S23C2S being the reduced 



with = J I sin (5 
Jarlskog invariant. The coefficients {a, 6, c} for P^^ can 
be obtained by replacing A with —A and J with —(5 in 
P^g. Hence, we obtain 

r. sin A A 

Aa ^ +8aJr0- — ^ — cos A, 

^ ^ sin A A . ^ 

Ab ^ -SaJ^e^ — smA, 



A 



Ac ^ 45^38^36+6 



(8) 



with 6± =sin[(A-l)A]/(A-l)±sin[(A-hl)A]/(A + l). 
Note that Eqs. ^ and ([8| are valid as long as a A <C 1, 



i.e., when the distance L and energy E are far away from 
the region where the low-frequency oscillation becomes 
dominant. This condition is satisfied in all the ongoing 
and upcoming long-baseline experiments, however, it is 
violated for atmospheric neutrino experiments. When 
low-frequency oscillations come into play, one can expand 
the probabilities in terms of 8^3, which are exact with 
respect to a, as is performed in Ref. [9 . 

Now, we apply the approximate formulas to the mea- 
sures of leptonic CP violation in Eqs. ([2|-([6| and explore 
their main features. 



A. CP Asymmetry A^^ (S) 

From Eqs. ^ and ([8|, one can obtain the approximate 
result for A^f[S). For neutrino energies E > 2 GeV, 
i.e., in the high-energy approximation, we can safely ig- 
nore Amli. In this two-fiavor limit, P^g is just given 
by the (5- independent term c. Note that the (5-dependent 
terms arise from the interference between the two mass- 
squared differences, and thus are suppressed by a factor 
a/5i3, indicating that the main structure of P^g in the 
L-E plane is determined by high-frequency and paramet- 
ric resonances [7 . Given the neutrino mass hierarchy, 
the resonance existing in the neutrino channel should be 
absent in the antineutrino channel, and vice versa. For 
lower energies, one has to consider three-ffavor oscilla- 
tions and analyze the resonance structure due to Am2i. 
For a systematic study of neutrino oscillograms of the 
Earth and the theoretical interpretation of their reso- 
nance structures, see Refs. 0^,. 

In Fig. 1, we have calculated A^f{5) by using the ex- 
act probabilities in the three-ffavor framework and the 
PREM model of the Earth matter density [17 . The base- 
line can be calculated via L = 2Rcos h with R :^ 6370 km 
being the Earth radius and h the nadir angle. The abrupt 
transition around h ~ 33° in Fig. 1 is caused by the 
change from the core- and mantle-crossing trajectories. 

The main structure in Fig. 1 is determined by mat- 
ter effects. In particular, the plots for 5 = indicate 
pure fake CP violation. In high-energy region, we have 
^^e" ^ 45^3S236^6_ to leading order. In the neutrino 
channel, the high-frequency resonance exists for the nor- 
mal mass hierarchy (NH). In other words, sin(74 — 1)A is 
resonantly enhanced in the NH case, while sin(74 + 1)A 
in the inverted mass hierarchy (IH). Therefore, from the 
NH case to the IH case, 6_ changes sign while 6^ not, 
which can be used to explain the sign-ffipping difference 
between the NH and IH cases. In Fig. 1, the numerical 
results are presented for the NH and IH cases in the up- 
per and lower rows, respectively. The sign-ffipping differ- 
ence can be clearly observed in the high-energy region. 
In the sub-GeV energy region, the low-frequency reso- 
nances and parametric enhancement comprise the dom- 
inant structure [8 . The peak and valley features come 
from the mismatch between neutrino and antineutrino 
probabilities. Hence, the results of A^f{S) in this re- 
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FIG. 1: Numerical results of A^^ {^) for 6 = and 7r/2 in the case of normal neutrino mass hierarchy (upper row) and inverted 
neutrino mass hierarchy (lower row), where the best-fit values O12 = 34°, Oi^ = 9°, O23 = 40°, Arn^i = 7.5 x 10~^ eV^, and 
lAmiil = 2.5 X 10-^ eV^ have been used p^. 



gion are essentially independent of the different neutrino 
mass hierarchies, and tiny differences may arise from the 
high-order corrections of Si^. The difference between the 
two plots in each row of Fig. 1 is more evident in the low- 
energy region, where the condition for the high-frequency 
resonance is not fulfilled. 



B. Working Observables AP^^ (6) and AP^^ 

Now, we turn to the working observables AP^^{S) and 
AP^ based only on P^^. Using Eqs. ^ and Q, we 
obtain the approximate formula 



Ap; 



CP 



-SaJ 



sinAA sm{A - 1)A 



A 



tan - cos A 
2 



A-1 
sin A 



(9) 



Note that the Jarlskog invariant J appearing here is 
just a simple notation and does not mean any difference 
between neutrino and antineutrino oscillations. On the 
other hand, from Eqs. Q and ([7|), we find 



AP" 



WaJ 



sinAA sin(A- 1)A 



A 



A 



(10) 



Since both AP^^ (S) and AP^ have already been stud- 
ied in great detail in Refs. [71 18], we will not repeat the 



analysis here. However, in Fig. 2, we have performed nu- 
merical calculations of AP^ for neutrinos and AP^ for 
antineutrinos in both NH and IH cases, to illustrate the 
main features. 

It is straightforward to understand the similarity be- 
tween the case of neutrinos with NH and that of antineu- 
trinos with IH. Note that AP^ for IH can be obtained 
by changing A —A, A —A, and a —a. On the 
other hand, the formula of AP^ for NH can be derived 
by replacing A —A. One immediately observes that 
these two formulas are identical. This feature is only 
present in the high-energy region, as shown in Fig. 2, 
where the similarity between the case of neutrinos with 
IH and that of antineutrinos with NH is also evident. 

There is a resonance deep in the core region for NH. 
But this resonance {h - 22° and 2 GeV < £; < 3 GeV) 
disappears for IH. The oscillatory structures in the low- 
energy region are different for neutrinos and antineutri- 
nos, but similar for NH and IH. Moreover, there are three 
vertical lines, i.e., the "solar magic lines" correspond- 
ing to APj^ = [8]. Solving sin A A = 0, one obtains 
AA = nTT, or L = 2n7r/V, with n being a positive inte- 
ger. On the other hand, the equation sin(A — 1)A = de- 
termines the "atmospheric magic lines". But these con- 
ditions for magic lines depend on the factorization ap- 
proximation (i.e., a ^ and s^^ -> 0), and the general 
forms can be found in Ref. [8]. In the low-energy region, 
the "solar magic lines" are no longer energy-independent, 
and they coincide with the low- frequency oscillation dips. 



5 




Working Observables AA^^ (6) and AA^^ 



Finally, we come to the new pair of working observ- 
ables AA^f{S) and AAJ^g. The motivation to introduce 
these observables is two- fold. First, any quantity measur- 
ing leptonic CP violation in neutrino oscillations should 
reflect the intrinsic difference between neutrino and an- 
tineutrino probabilities. Second, it should be helpful for 
an optimal experimental setup to practically measure 6. 
Therefore, we derive the working observables from A^f . 
Combining Eq. ([5| with Eq. ([8|, we obtain 



we flnd 



A .rp ^ sin A A 



6. 



tan - cos A 
2 



e+ sin A 



which reduces to the vacuum result A^f when A 



(11) 
-> 0. 

As AA^f is proportional to J7, it vanishes for ^ = and 
the CP-violating effects induced by matter effects have 
been removed. In Fig. 3, we show the numerical results 
of AA^f for (5 = 7r/2 in both NH and IH. By deflni- 
tion, they are just the differences between the two plots 
in each row of Fig. 1. Apart from a resonance region in 
the deep core, sizable values of AA^f are lying in the 
energy region below 1 GeV. Note that the results for NH 
and those for IH in the high-energy region are nearly in- 
distinguishable, which can be understood by noting that 
Eq. (11) is invariant under the transformations a —a, 
A -A, and A -A. 

To examine the sensitivity to we consider the range 
of A^^{6) by varying 6 in [0, 27r). From Eqs. (p| and (pi. 



:sin A A 



AA"^^ = IGaJ^y (©-CosA)2 + (6+ sin A)' 

(12) 

Similar to AP^, the domain structure of AA^f (S) and 
AAJ^g can be understood through the "solar" and "at- 
mospheric" magic lines, and the interference phase con- 
dition, as suggested in Ref. [8 . In the upper (lower) row 
of Fig. 4, we give the numerical results of AAJ^^ in NH 
(IH). It is now evident that the most signiflcant value 
appears in the area of relatively short baselines and low 
energies, which are quite relevant for the long-baseline 
experiments. The zoom-in plots of this region are shown 
in the right column of Fig. 4, where the ongoing and 
upcoming neutrino experiments have also been indicated 
by solid diamonds and the energy ranges are represented 
by the peak energies plus error bars. It is worthwhile to 
mention that AAJ^g in Eq. (12) is unchanged when we 
switch from NH to IH through a^— a, A^— A, and 
A —A. Hence, AAJ^^ is insensitive to the neutrino 
mass hierarchy. 



D. Optimal Experimental Setup 

Now, we discuss the optimal experimental setup to 
probe S. The ongoing long-baseline experiments T2K 
{E = 0.72 ± 0.27 GeV, L = 295 km) and NOuA 
{E = 2.02 ± 0.43 GeV, L = 810 km), together with 
the proposed ones LBNE {E = 3.55 ± 1.38 GeV, L = 



6 




FIG. 3: Numerical results of AA^f{5) for 5 = 7r/2 in the cases of normal neutrino mass hierarchy (left) and inverted neutrino 
mass hierarchy (right), where the input neutrino parameters are the same as in Fig. 1. 




h [1 h i°] 

FIG. 4: Numerical results of AA™e in the case of normal neutrino mass hierarchy (upper row) and inverted neutrino mass 
hierarchy (lower row), where the input neutrino parameters are the same as in Fig. 1. In the right plots, we have zoomed in 
the parameter region that is relevant for the long-baseline neutrino oscillation experiments T2K (green, at h — 88.7°), NOz/A 
(cyan, at /i = 86.4°), LBNE (red, at /i = 84.1°), LAGUNA-LBNO (blue, h = 79.7°), and ESS (black, at /i = 87.8°). 



1300 km) and LAGUNA-LBNO^ {E = 5.05 ± 1.65 GeV, 
L = 2288 km), will be considered for illustration. Ex- 
cept for T2K, all experiments are intended to equally 
operate both in the neutrino and antineutrino channels. 
Therefore, one can construct the CP asymmetry A^f 
by measuring the neutrino and antineutrino probabili- 
ties. In the right column of Fig. 4, we observe that these 
experimental setups are lying on the "first band" (i.e.. 



At the moment, the fate of this project is unclear. 



AAJ^g ~ 10 %), counted from top-right to bottom- left. 

To improve the experimental sensitivity, one can lower 
the neutrino beam energy and locate the experiment on 
the "second band" (i.e., AAJ^g ~ 15 %). However, this 
observation is based on the probability level and we have 
to notice the energy dependence of the neutrino flux and 
the cross section. In addition, the detection efficiency and 
the background should be taken into account. Therefore, 
a detailed simulation has to be performed in order to 
make a final conclusion. Recently, it has been proposed 
that the ESS proton beam can be adjusted to produce 
an intense neutrino beam of energy around 0.4 GeV [ 18] . 
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The simulation results indicate that with 8 years of data 
taking with an antineutrino beam and 2 years with a 
neutrino beam up to 73 % of the whole range of 5 could 
be covered at 3a level at the optimal baseline of around 
500 km [TSl. In the right column of Fig. 4, the ESS 
proposal is shown as £^ = 0.4 ±0.2 GeV and L = 500 km. 
Such an experimental setup happens to be on the "second 
band" , as we suggest. For preliminary performance of the 
different experimental setups, see Refs. p!8H23] 

It is worthwhile to mention that precision measure- 
ments of atmospheric neutrinos in the PINGU detector 
at the IceCube experiment on the South Pole may have 
good sensitivity to S [24 . A study of the PINGU de- 
tector and accelerator-based neutrino beams also exists 
[25 . On the other hand, there is no doubt that the future 
neutrino factories are the best place to measure S with 
high statistical significance [26H29] . 

IV. SUMMARY 

We have made a complete survey of measures of lep- 
tonic CP violation in neutrino oscillation experiments. 
Two new working observables AA^^{S) = A^^{S) — 
ACP(O) and AA-^ = me^[A'^^{d)]-min[A'^^{d)] , where 

^^ei^) = Pfiei^) ~ ^fiei^) asymmetry of oscilla- 

tion probabilities, are suggested to describe the intrinsic 
leptonic CP violation. Both analytical and numerical cal- 



culations are performed to illustrate their main features. 
The band structure of AA^^ in the baseline-energy plane 
can be implemented to optimize the experimental setup. 
Furthermore, we have found that the current and future 
long-baseline neutrino oscillation experiments are located 
on the first band. On probability level, we observe that 
the decrease of the neutrino beam energy in a proper way 
could improve the experimental sensitivity. 

Although the final verdict on the optimal experimen- 
tal setup requires a more sophisticated simulation, we 
expect our analysis to be helpful in understanding the 
leptonic CP violation and useful in probing the leptonic 
CP-violating phase in future long-baseline experiments. 
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